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Abstract 

We update the upper bound on the lightest CP even Higgs mass in the NMSSM, 
which is given as a function of tan/3 and A. We include the available one and two loop 
corrections to the NMSSM Higgs masses, and constraints from the absence of Landau 
singularities below the GUT scale as well as from the stability of the NMSSM Higgs 
potential. For mt op varying between 171.4 and 178 GeV, squark masses of 1 TeV and 
maximal mixing the upper bound is assumed near tan/3 ~ 2 and varies between 139.9 
and 141.4 GeV. 
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1 Introduction 



Supersymmetric extensions of the standard model predict quite generally at least one rela- 
tively light Higgs boson. Hence, as soon as results from future collider experiments provide us 
with informations on the mass of at least one Higgs boson, we will be able to put constraints 
on possible supersymmetric extensions of the standard model. 

To this end we need to know, as accurately as possible, how the Higgs boson masses 
depend on the nature of the supersymmetric extensions of the standard model, and on the 
parameters of these models. (We hope, of course, to get independent informations on these 
parameters from direct sparticle detections in the future.) 

In the MSSM, corresponding calculations have been pushed to a fairly high accuracy, 
including many two-loop corrections. Recent reviews on the lightest Higgs boson mass in 
the MSSM can be found in refs. [1-3]. 

In the present paper we discuss the simplest version of the NMSSM [4] with a scale 
invariant superpotential 

W = \SH u H d + ^S 3 + ... , (1.1) 

which is the only supersymmetric extension of the Standard Model where the weak scale 
originates from the soft susy breaking scale only, i.e. where no supersymmetric dimensionful 
parameters as are present in the superpotential. 

It is well known [4] that the lightest Higgs boson in the NMSSM can be heavier than 
the one of the MSSM due to additional terms in the tree level Higgs potential proportional 
to A 2 ; the additional contribution is 

\2 

Am 2 = — Mf sin 2 2/3 . (1.2) 

If one requires the absence of a Landau singularity for A below the GUT scale, A is bounded 
by ~ 0.7 from above [4], leading still to an upper bound on the mass of the lightest Higgs 
boson that is, however, larger than in the MSSM. Thus, future measurements of the Higgs 
boson mass could serve to distinguish these two models, provided that we know the difference 
between the upper bound on the mass of the lightest Higgs boson in the different models. 

At present, the radiative corrections in the NMSSM have not been computed to quite 
the same accuracy as in the MSSM. Of course, radiative corrections in the NMSSM that 
are proportional to the quark/lepton Yukawa couplings and the gauge couplings only are 
the same as in the MSSM, but there are many additional contributions involving the new 
Yukawa couplings A and k in the superpotential in eq. (1.1), and the associated soft trilinear 
couplings A\ and A K . 
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The one loop corrections in the NMSSM induced by t and b quark/squark loops have 
been computed already some time ago [5], and the dominant two loop corrections (~ h\ 
and ~ hja s ), that are the same as in the MSSM, have been included in an analysis of the 
NMSSM Higgs sector in ref. [6, 7]. 

The leading logarithmic one loop corrections to the lightest Higgs mass in the NMSSM 
proportional to the electroweak gauge couplings g or NMSSM specific Yukawa couplings A, 
K (~ Q^iQ 2 ^ 2 ^ 2 ^ 2 ) have been computed only recently [8]. They are included in the 

latest version of the code NMHDECAY [9-11], where the NMSSM Higgs masses, couplings 
and branching ratios are computed as functions of the parameters in the Lagrangian of the 
model. 

This code checks also the absence of a Landau singularity for A below the GUT scale 
using the two loop renormalization group equations, and susy threshold effects around the 
susy scale. This procedure is numerically relevant, since Am\ depends on A, and the upper 
bound on A depends on tan/5 (via the top Yukawa coupling h t ) and k. 

It is the purpose of the present paper to review the upper bound on the lightest Higgs 
boson mass in the NMSSM, using the up-to-date knowledge of the corresponding radiative 
corrections. 

Instead of investigating far-fetched regions in parameter space that serve to obtain very 
conservative bounds, we proceed as follows: For the soft terms that are relevant for the 
sparticle spectrum, we chose universal squark and slepton masses of 1 TeV, and trilinear 
couplings of 2.5 TeV (that practically maximize the one loop radiative corrections to the 
Higgs boson mass). For the gaugino masses we take Mi = 150 GeV, M 2 — 300 GeV and 
M3 = 1 TeV in rough agreement with universal gaugino masses at the GUT scale. For the 
top quark pole mass we present results both for m top = 171.4 GeV (the latest central value 
obtained by the Tevatron Electroweak Working Group [12]) and a very conservative upper 
limit of m top = 178 GeV. The NMSSM specific Yukawa couplings and trilinear soft terms 
A, k, A\ and A K as well as the effective \x parameter (// = As in the NMSSM) are chosen 
such that the lightest Higgs boson mass is maximized, without violating constraints from 
perturbativity of the Yukawa couplings at the GUT scale, nor phenomenological constraints 
on CP odd or charged Higgs masses and couplings. To this end a numerical analysis is 
required, that is performed using the updated version of NMHDECAY [11]. The upper bound 
on the lightest Higgs mass is then given as a function of tan/5 and A. 

For the same choice of the above soft terms, we present the upper bound on the lightest 
Higgs mass in the MSSM limit A — > as obtained with NMHDECAY. This result can be 
compared to values obtained from analytical or numerical analyses in the MSSM, that include 
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radiative corrections that are still absent in NMH DECAY: these are notably electroweak one- 
loop corrections ~ g 4 beyond the LLA, and non-dominant two-loop corrections (involving 
less than two powers of large logarithms) ~ and ~ h\ a s beyond the ones that follow 
from an RG improvement of the one loop corrections [13] (which are included). 

For soft terms as above, rrit op = 178 GeV, Ma = fi = 1 TeV and tan/5 = 10 SuSpect 
gives rrih ~ 128.5 GeV (taken from ref. [2]), FeynHiggs rrih ~ 134 GeV (taken from ref. [3]), 
and NMHDECAY rrih ~ 128.6 GeV. This allows to estimate the uncertainties on rrih due to 
the radiative corrections not included in NMHDECAY, following the discussions in [2,3] that 
we will not repeat here. 

The striking effect in the NMSSM is that the maximal value of rrih is not assumed for large 
tan/5 as in the MSSM, but at low tan/5 ~ 2 due to the tree level term noted above. There 
we obtain rrih ~ 139.9 GeV for m top = 171.4 GeV, and rrih ~ 141.4 GeV for mt op = 178 GeV 
(for the same other parameters as above). For larger values of tan/5 the upper bound on rrih 
decreases in the NMSSM. For tan/5 > 10 it hardly exceeds the MSSM value given above, 
since the effect of the tree level term becomes small. For small tan/5 < 2 the absence of a 
Landau singularity below Mqvt restricts A more strongly from above, due to the large top 
Yukawa coupling ht. This implies that present lower limits on rrih from LEP still lead to a 
lower bound on tan/5 of ~ 1.3 in the NMSSM. 

It must be noted that in particular regions of the parameter space of the NMSSM the 
upper bound on rrih discussed here can be misleading: 

In principle, a singlet-like CP even Higgs boson can be lighter than the lightest doublet- 
like CP even Higgs boson (with non-vanishing couplings to the Z boson) in the NMSSM. 
Strictly speaking, the upper bound on the lightest CP even Higgs boson discussed here is 
then still valid. 

However, a singlet-like CP even Higgs boson would have been practically undetectable 
at LEP due to its vanishing coupling to the Z boson. Fortunately, if the lightest CP even 
Higgs boson is a pure singlet in the NMSSM, the upper bound on rrih discussed here applies 
then to the lightest doublet-like CP even Higgs boson. On the other hand, if the lightest CP 
even Higgs boson is only approximately a singlet, the lightest doublet-like CP even Higgs 
boson can be heavier than the upper bound on rrih discussed here. 

A similar reasonning applies to the situation where the doublet-like CP even (SM like) 
Higgs boson decays into singlet like (mostly two CP odd) scalars [14]. Then, the detection of 
the SM like Higgs boson can be very challenging, even if its mass satisfies the upper bounds 
discussed here. Hence, although the upper bound on rrih presented here is always valid, it 
may refer to a state that is difficult to detect. 
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2 The upper bound on the lightest Higgs Boson mass 
in the NMSSM 



In order to find the regions in the parameter space of the NMSSM that maximize the upper 
bound on the lightest CP even Higgs Boson mass, it is helpful to take a look at the CP 
even Higgs mass matrix at tree level. In the basis (H u , H d , S) and using the minimization 
equations in order to eliminate the soft masses squared, it reads: 



M 
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9 2 K + ^ (A\ + v) (2A 2 - g 2 )h u h d - fi(A\ + v) 2\h u pi - \h d (A x + 2v) 

h 

g 2 h 2 d + fi-^(A x + u) 2\h df i - Xh u (A x + 2v) 

d h h 

1 2 A ^ul^d 



\ 2 A X -^ +v(A K + 4v) 
/i 

(2.1) 

where v = ks. To a good approximation, the 2x2 doublet subsector is diagonalized by 
the angle (5 which gives the desired light eigenstate h and a heavy eigenstate H with a 
mass m H ~ nriA close to the MSSM-like CP odd state (the larger tua, the better this 
approximation). In the NMSSM, one can define m\ as the diagonal doublet term in the CP 
odd 2x2 mass matrix after the Goldstone mode has been dropped. At tree level, it has the 
same expression as in the MSSM: 

m A = , wi th B = A A + v . (2.2) 

sin 2p 

In the CP even sector of the NMSSM this is not the end of the story, however: the light 
eigenstate h of the 2x2 doublet subsector still mixes with the singlet state S, which is 
heavier than h by assumption. In order to maximize mh, this mixing has to vanish: 

A [2n - (A x + 2v) sin 2(5] ~ . (2.3) 

This requires either A — > (which minimizes the NMSSM specific tree level contribution to 
m h of eq. (1.2)) or 

Ax ^ - 21/ . (2.4) 

sin 2p 

On A K we get the following constraints: -M| 33 in eq. (2.1) must at least be positive, 
which requires essentially (since the first term is typically relatively small) A K v > — 4z/ 2 . 
The CP odd mass matrix element in the singlet sector, given by 

M 2 Pi33 = A\Kh u h d + \ 2 A\^—^ - 3uA K , (2.5) 



5 



must also be positive. Typically the last term in (2.5) dominates, hence we get an allowed 
window 

-4z/ 2 < A R v < . (2.6) 

Next, in order to maximize the NMSSM specific tree level contribution to m h of eq. (1.2), 
A has to be as large as possible; we require, however, the absence of a Landau singularity for 
all Yukawa couplings A, k, h t below the GUT scale, which leads to the following constraints: 

First, given the corresponding RG equations [4], this implies small values for k. The limit 
k — ► while A remains finite is disallowed, however, both from the stability of the potential 
and the fact that the allowed window of eq. (2.6) vanishes in this limit. (Of course, stability 
of the potential and positivity of all masses squared are related issues.) 

Second, for small tan/3 the top quark Yukawa coupling becomes large, and can run into 
a Landau singularity below Mqut, or induce a Landau singularity below Mqut for A. The 
value of tan/3 that allows for maximal values of A (and maximizes the tree level contribution 
to m h of eq (1.2)) is around 2. 

In this region of tan/3, a larger value for the top quark pole mass does hardly increase 
the upper bound on m^: at fixed tan/3, larger m top implies a larger top Yukawa coupling 
h t , which implies a somewhat lower allowed value for A. Consequently a variation of the 
top quark pole mass between 171.4 and 178 GeV (which increases m h by ~ 4.8 GeV for 
tan/? ~ 10), increases the maximal allowed value for m h in the NMSSM by only ~ 1.5 GeV 
for tan/3 ~ 2. 

For large values of tan/3, it is obvious from eq. (2.2) that m A tends to be very large, 
unless B is small, (pi cannot be smaller than ~ 100 GeV due to the lower bound on chargino 
masses from LEP). Very large values of uia are unnatural, since they require supersymmetry 
breaking Higgs masses of the same order of magnitude, which aggravate the fine tuning 
problem - a situation which we want to avoid. In the MSSM, one can always chose B small 
enough to keep m A reasonable even at large tan/3. In the NMSSM this is also possible, 
provided that v ~ — A\. However, one has also to minimize the doublet-singlet mixing of 
eq. (2.3) in order to to maximize m h . If A is not very small, eq. (2.4) together with v ~ — A x 
implies z/^>/i, which is equivalent to k^$> A. Large values of k leading to a Landau singularity 
below the GUT scale, this is excluded. Thus, the only way of minimizing the doublet-singlet 
mixing while keeping uia constant at large tan/3 is to assume A — > 0, which means that the 
bound on m h is the same as in the MSSM. 

(In general, for large values of tan/3 the LEP constraints on m h imply either that m A and 
\A\\ assume very large values > 1 TeV, or A < 0.2.) 
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All these considerations make it clear that a realistic upper limit on rrih in the NMSSM 
requires numerical methods; analytic approaches can be misleading (and can allow for larger 
values of rrih). 

Our results below are obtained with NMHDECAY [11]. The precision of the included 
radiative corrections to the lightest CP even Higgs mass has already been discussed in the 
introduction and is given in [9] and [10]. 

As discussed in the introduction, we take universal squark and slepton masses of 1 TeV, 
and trilinear squark/slepton couplings of 2.5 TeV (near maximal mixing). For the gaugino 
masses we take M x = 150 GeV, M 2 = 300 GeV and M 3 = 1 TeV. We scan over the NMSSM 
specific Yukawa couplings and trilinear soft terms A, re, A\ and A K as well as the effective fx 
parameter, and we obtain the regions in the NMSSM parameter space that maximize rrih in 
agreement with the considerations above. 

In fig. 1 we show our results for the the upper bound on rrih for 1 < tan/? < 10. The 
thick full line corresponds to m top = 178 GeV, the thin full line to m top = 171.4 GeV, both 
without imposing constraints on rriA- 

With the above soft terms, the upper bound on rrih in the NMSSM is 141.4 GeV for 
m top = 178 GeV. It is reached for tan/? ~ 2.2, A ~ .677, re ~ .068, fx ~ 545 GeV, A\ ~ 
1365 GeV, and A K ~ 10 GeV (strictly speaking a certain range of values for re, fx, Ax and 
A K gives the same result for rrih for these values of tan/? and A). For m top = 171 .4 GeV, 
the upper bound on rrih is 139.9 GeV and is obtained for tan/3 ~ 2, A ~ .703, re ~ .049, 
fx ~ 534 GeV, A x ~ 1287 GeV and A K ~ 10 GeV. 

For tan/? = 10 we get 133.6 GeV in the NMSSM for m top = 178 GeV (resp. 128.8 GeV for 
iTitop — 171 .4 GeV), which remains nearly constant for larger values of tan/? (a slight increase 
of the contributions from the radiative corrections is compensated by a slight decrease of the 
tree level term of eq. (1.2)). 

In the same fig. 1, we show the upper bound on rrih m the MSSM limit A — > as obtained 
with NMHDECAY as a thick dashed line for m top = 178 GeV, and as a thin dashed line for 
mtop — 171.4 GeV (taking rriA = 1 TeV). In this limit, the upper bound on rrih reaches 
129.7 GeV for m top = 178 GeV (resp. 124.4 GeV for m top = 171.4 GeV) at tan/? = 10, and 
increases by another 1 GeV for very large tan/? = 50. 

As noted above, large values of tan/? imply large values for rriA in the NMSSM, if A is 
kept fixed. Indeed, along the full lines of fig. 1 the value of rriA increases with tan/? up to 
several TeV. The consequence of fixing rriA < 1 TeV is that the maximally allowed value 
of A decreases with tan/?. The corresponding effect on the upper bound of rrih is shown as 
a thick dotted line for m top = 178 GeV, and a thin dotted line for m top = 171 A GeV in 
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Figure 1: Upper bound on the lightest Higgs mass in the NMSSM for m top = 178 GeV 
(thick full line: arbitrary, thick dotted line: = 1 TeV) and m top = 171.4 GeV 

(thin full line: arbitrary, thick dotted line: = 1 TeV) and in the MSSM (with 
rfiA = 1 TeV) for m top = 178 GeV (thick dashed line) and m tov = 171.4 GeV (thin dashed 
line) as obtained with NMH DECAY as a function of tan/3. Squark and gluino masses are 
1 TeV and A top = 2.5 TeV. 

fig. I. Now we get an upper bound of 130.1 GeV for m top = 178 GeV (resp. 124.7 GeV for 
m to P = 171.4 GeV) at tan/3 = 10. For larger values of tan/3, the upper bound on m h remains 
essentially the same as in the MSSM. 

Hence, our main result is that the upper bound on m h is ~ 12 GeV (for m top = 178 GeV) 
or ~ 16 GeV (for m top = 171.4 GeV) larger in the NMSSM as compared to the MSSM, and 
is obtained for small tan/3. For very large tan/3, the difference between the upper bound on 
m h in NMSSM and in the MSSM vanishes, if tua is assumed to remain smaller than a few 
TeV. 

Let us compare this bound on rrih to earlier work: it is about 6 GeV larger than the 
one obtained from fig. 4 in ref. [7] (for the corresponding values for m top ). Also the value of 
tan/3, where this bound is reached, is now smaller (~ 2 compared to ~ 3 in ref. [7]). These 
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Figure 2: Upper bound on the lightest Higgs mass in the NMSSM for m top = 178 GeV, 
tan/3 = 2.2, electroweak/Yukawa corrections included (thick full line) and omitted (thick 
dotted line), and m top = 171.4 GeV, tan/3 = 2, electroweak/Yukawa corrections included 
(thin full line) and omitted (thin dotted line). Squark and gluino masses and A top are as in 
fig. 1. 

differences are due to the improved treatment of radiative corrections in NMH DECAY which 
concerns both the two loop corrections ~ h\ and h\a a (which are now RG-improved) , and the 
inclusion of one loop corrections (in the LLA, keeping terms ~ ln(M| /M§)) proportional 
to the electroweak gauge couplings and NMSSM specific Yukawa couplings A and k. The 
effect of the first improvement is a considerable increase in rrih, whereas the effect of the 
electroweak/Yukawa corrections is a slight decrease of m h by up to ~ 2 GeV. 

In order to clarify the latter effect and, simultaneously, the general effect of the NMSSM 
specific Yukawa couplings at low tan/3, we show in fig. 2 the upper bound on rrih as a function 
of A at fixed tan/3. Here the thick full line corresponds to m top = 178 GeV, tan/3 = 2.2 and 
electroweak/Yukawa corrections included, whereas the thick dotted line would be the result 
with these corrections omitted. The thin full line corresponds to m top = 111 A GeV, tan/3 = 2 
and electroweak/Yukawa corrections included, whereas the thin dotted line would be the 
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result without these corrections. One sees the decrease in rrih due to the electroweak/ Yukawa 
corrections, which increases the lower bound on A for small values of tan/5 and m top due to 
the LEP bound on nih- 

As final remark we repeat, as noted at the end of the introduction, that the mass of the 
lightest detectable Higgs boson could be larger in the NMSSM than the upper bounds given 
here; in order to interpret future data in the context of the NMSSM, constraints (or positive 
results) must be available in the plane Higgs mass versus Higgs couplings in order to be 
sensitive to a possible singlet /doublet mixing. 
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